Magnetic incommensurability in p-type cuprate perovskites 



A. Sherman 

Institute of Physics, University of Tartu, Riia 142, 51014 Tartu, Estonia 

(Dated: August 23, 2011) 

Abstract 

For the superconducting phase with a d-wave order parameter and zero temperature the mag- 
netic susceptibility of the t-J model is calculated using the Mori projection operator technique. 
Conditions for the appearance of an incommensurate magnetic response below the resonance fre- 
quency are identified. A fast decay of the tails of the hole coherent peaks and a weak intensity 
of the hole incoherent continuum near the Fermi level are enough to produce an incommensurate 
response using different hole dispersions established for p-type cuprates, in which such response 
was observed. In this case, the nesting of the itinerant-electron theory or the charge modulation of 
the stripe theory is unnecessary for the incommensurability. The theory reproduces the hourglass 
dispersion of the susceptibility maxima with their location in the momentum space similar to that 
observed experimentally. The upper branch of the dispersion stems from the excitations of localized 
spins, while the lower one is due to the incommensurate maxima of their damping. The narrow 
and intensive resonance peak arises if the frequency of these excitations at the antiferromagnetic 
momentum lies below the edge of the two-fermion continuum; otherwise the maximum is broad 
and less intensive. 

PACS numbers: 74.72.-h, 75.10.-b, 75.40.Gb 
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I. INTRODUCTION 



One of the most interesting features of the inelastic neutron scattering in lanthanum 
cuprates is that for hole concentrations x > 0.04, low temperatures and small energy 
transfers the scattering intensity is peaked at the incommensurate momenta f|, |±«J, 
(| ± k, |) in the reciprocal lattice units 27r/a with the lattice period a.—^ For x < 0.12 
the incommensurability parameter k is approximately equal to x and saturates for larger 
concentrations.- The incommensurability was observed both below and above the temper- 
ature of the superconducting transition T c .->^ The analogous low-frequency incommensu- 
rability was observed also in YBa 2 Cu 3 7 _, / and Bi 2 Sr2CaCu 2 8 ." This gives ground to 
suppose that the incommensurability is a common feature of p-type cuprate perovskites. 
As the transfer frequency increases, the incommensurability parameter k decreases and at 
the frequency u r ~ 25 — 40 meV the magnetic response becomes commensurate, being 
peaked at the antiferromagnetic wave vector Q = (|,|)--^ The value of the u r depends 
on the hole concentration. For even larger frequencies in some crystals the magnetic re- 
sponse becomes again incommensurate, with the incommensurability parameter growing 
with frequency.- 1 ^ In some experiments susceptibility maxima were found at momenta 
(I ±«, | ±k), f| ± k, | =F ft).-~— In other experiments the maxima form a circle around 
Q or even merge together into a broad commensurate maximum.— ~— Thus, the dispersion 
of the maxima in the susceptibility resembles a hourglass with the neck at the momentum Q 
and the frequency cj r ."^ In YBa2Cus07_ y , in the region of the neck the susceptibility peak 
is much more intensive and narrower than the maxima for smaller and larger frequencies. 
This peculiarity of the spin-excitation spectrum was called the resonance peak.— 

The most frequently cited theories for the magnetic incommensurability are based on 
the picture of itinerant electrons with the magnetic susceptibility calculated in the random 
phase approximation^— and on the stripe picture.— >22rM In the former approach the in- 
commensurability below bj r is connected with a nesting in the constant energy surfaces of 
the carrier dispersion. The approach allows one to reproduce the observed hourglass disper- 
sion of the susceptibility maxima. However, carrier dispersions derived from photoemission 
experiments frequently have no nesting,— and some fitting of their parameters is neces- 
sary to obtain the incommensurability.- 1 Such sensitivity of the results to fine details of the 
hole dispersion casts some doubt upon this approach, since in different p-type cuprates, in 
which the incommensurability is observed, the dispersions diverge considerably and vary 
with doping. Another theory^ based on the itinerant-electron approximation relates the 
incommensurability to a singularity in the noninteracting susceptibility, which arises on the 
assumption of vanishingly small damping of fermion states.— 1 ^ 

The second approach is based on the notion of charge stripes - a one-dimensional pe- 
riodic variation of the charge density in a Cu-0 plane. This theory also reproduces the 
observed hourglass dispersion of the susceptibility maxima with their proper location in the 
Brillouin zone, at least for x — |. It should be noted that static charge stripes are ob- 
served only in crystals with the low-temperature tetragonal (LTT) or the less-orthorhombic 
phases (L^-y-^NdySr^CuC^ and La 2 _ a; Ba a ;Cu04), and such stripes are not detected in the 
low-temperature orthorhombic (LTO) phase of cuprates, e.g., in La2_a.Srj.CuO4.— ^ At the 
same time the magnetic incommensurability is similar in these phases. To resolve this con- 
tradiction the supposition of stripe fluctuations in the LTO phase with the same symmetry as 
in the LTT phase was proposed.— However, there are strong grounds to believe that stripes 
are connected with deformations of Cu-0 planes due to certain tilts of Cu06 octahedra.— ~— 
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In this case in the LTO phase dominant fluctuations will be oriented along diagonals of the 
planes.— With such fluctuations the susceptibility maxima below u r will be located along 
diagonals of the Brillouin zone, in contrast to experimental observations. Besides, at fre- 
quencies u > u r the stripe theory predicts maxima of the magnetic susceptibility, which are 
sharply peaked in the momentum space at (~ ± n, ~ ± k), (| ± k, ~ =f k). In experiments 
the susceptibility looks more isotropic in this region.— 1 ^ 



Strong electron correlations in Cu-0 planes were taken into account in Refs. 34.35 
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calculating the magnetic susceptibility of the two-dimensional (2D) t-J model. To carry out 
calculations with the Hubbard operators entering into the Hamiltonian the Mori projection 



operator technique^ was used. In Ref. |34j, attention was drawn to the fact that for small 
frequencies the incommensurate response points to a dip in the spin-excitation damping at 
Q. The dip was related to the nesting of small pockets around (±|,±^J, which form the 
Fermi surface of the t-J model at small x. The nesting is a consequence of the short-range 
antiferromagnetic order with a large correlation length, which is established for small hole 
concentrations. Some similarity between equations for the susceptibility obtained in the Mori 
and itinerant-electron approaches was pointed out in Ref. |35|. In this work, the arguments 
of the itinerant-electron picture were used for the explanation of the incommensurability 
below and above u r . In the case of strong electron correlations the incommensurability 
above the resonance frequency was related to the dispersion of excitations of localized spins in 
Refs. l37|]38l . In Ref. |39| this part of the hourglass dispersion was also obtained in the approach 
based on the Mori projection operator technique. However, due to a poor representation 
for the operator, which determines the self-energy of spin excitations, the contribution of 
the decay into an electron-hole pair was lost. Together with it the ability to describe the 
low- frequency incommensurability was gone. 

In the present work, we calculate the magnetic susceptibility of the 2D t-J model using 
the Mori projection operator technique. We take into account all contributions to the 
polarization operator of spin excitations. These contributions are connected with a decay of 
a spin excitation into a fermion pair, a decay into the pair assisted by a hole or with a creation 
of another spin excitation, and a decay into three spin excitations. In the calculations, we 
take into account an incoherent hole-spin-excitation (HSE) continuum, which holds the 
major part of the quasiparticle weight of the hole spectral function in the considered case 
of strong correlations. For the coherent part of the spectral function several dispersions 
found in the literature were analyzed, including those designed to fit photoemission data. 
For the superconducting state with the <i-wave order parameter Ak we found that at certain 
conditions all considered dispersions produce an incommensurate magnetic response below 
the resonance frequency. This differs from the result obtained in the itinerant-electron 
picture where, as mentioned above, a special fitting of the hole dispersion was necessary 
to obtain the incommensurability. The cause of the difference is in the following. As it 
will be seen below, for moderate hole concentrations and frequencies less than 150 meV, 
there is a certain similarity in the description of the magnetic susceptibility in the Mori 
and itinerant-electron approaches. In the latter theory, (^-functions in the imaginary part of 
the noninteracting susceptibility are usually approximated by Lorentzians with an artificial 
broadening. 19,21 An inspection of this approximation shows that far tails of the Lorentzians 
strengthen the intensity at Q. In this situation, the Fermi surface has to contain nearly 
parallel regions for the intensity at the respective incommensurate wave vector be larger than 
at the commensurate Q. The matter concerns tails, which are separated from the Lorentzian 
maximum by frequencies of the order or larger than the superconducting gap A = max(Ak), 
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i.e. more than an order of magnitude larger than estimated widths of spectral peaks. Such 
long tails imply that decay widths of states near the Fermi level are independent of frequency. 
However, as known, the damping grows with distance from the Fermi level.— ^ Therefore, 
the tails of the peaks decay much faster than in the Lorentzian. Besides, the tails have to 
be cut off in the superconducting gap. To take these facts into account, we approximated 
the (^-functions, which ensure the energy conservation, by a II-shaped step function of a 
finite width. This approximation gives the incommensurate magnetic response below the 
resonance frequency for the mentioned dispersions provided that the intensity of the HSE 
continuum is weak in the ranges |Ak| < \oo\ < A. 

This approach allows us to describe the hourglass dispersion of the susceptibility maxima, 
in which the part above u r reflects the dispersion of spin excitations. The maxima below 
u r originate from incommensurate peaks of the spin-excitation damping. They are sharply 
peaked at the axes of the Brillouin zone, while above the resonance frequency the suscepti- 
bility is more isotropic. Notice, however, that for one of the considered hole dispersions the 
low-frequency maxima were located on the diagonals, and in some cases the high-frequency 
incommensurate peaks merge into a broad commensurate maximum. The evolution of the 
Fermi surface and the increase of the spin-excitation frequency uq with doping explains 
rather naturally the observed growth of the incommensurability parameter k with x. In this 
approach, a narrow and intensive resonance peak arises in the susceptibility when the fre- 
quency uq falls on the region of small damping below the edge of the two-fermion states at 
the antiferromagnetic wave vector. If uq exceeds this threshold, a broader and less intensive 
maximum appears. 

The main formulas used in the calculations are given in the next section. In more details 
their derivation is presented in the Appendix. The discussion of the obtained results is 
carried out in Sec. III. Section IV is devoted to the speculation about possible changes in 
the susceptibility in the pseudogap and phase separated states and to concluding remarks. 

II. MAIN FORMULAS 

The Hamiltonian of the 2D t-J model reads 

H = X>oLa P( , + ~£ Jiv + sjfV) - A*Z)*i°> (1) 

ll'cr U' 1 

where a\„ = |lc)(10| is the hole annihilation operator, \la) and |10) are singly occupied and 
empty states on a site 1 of a 2D square lattice, which models a Cu-0 plane of p-type cuprate 
perovskites, a = ±1, t\y and J\y are hopping and exchange constants, s\ = \ J2 a cr l^ cr )(^ cr l 
and sf = | lcr) (1, — a\ are components of the spin-| operator, the Hubbard operator Aj 00 = 
1 10) (10 1, and fi the chemical potential. 

Our aim is the calculation of Green's function of spin excitations D(ku)), which is con- 
nected with the magnetic susceptibility by the relation^ 

X (M = -4^|D(kw) (2) 

with the Bohr magneton fis- The function reads 

/oo 
e**DQd;)dt, D(kt) = -i9(t) { [s z k (t), £j > . (3) 
-oo 
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In this equation, the averaging denoted by the angle brackets and the time dependence of 
the operator are defined by Hamiltonian (OQ), s£ = iV -1 / 2 e~* kl sf , and A" is the number 
of lattice sites. 

In the considered case of strong electron correlations it is convenient to use the Mori 
projection operator formalism^ for calculating Green's function ([3]). The respective formulas 
are given in the Appendix. The obtained expression for the function reads 



D(kio) 



hi 



to 2 — coU(kuj) — to 



2 ■ 



In this equation, 



/i k = 4(tF\ + Jd)(7k-1) 
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where we took into account only the nearest neighbor hopping and exchange terms in Hamil- 
tonian (P), which are characterized by the constants t and J, respectively. The multiplier a 
is introduced to correct the decoupling procedure^ 2 - 1 ^ (see the Appendix), 



a k <x = N- 1 ' 2 Y,i e m a la , s£ = N" 1 / 2 J2i e" iakl sf and 7 k = \ [cos{k x ) + cos^)]. Here and 
below the lattice spacing a is set as the unit of length. 

In Eq. (jlj), the polarization operator is given by the expression 



n(kw) 



-K 1 



N 2 



£ [2/i 2 (kqq') + /|(kqq')] / / / du l( ku 2 dw 3 



qq 



x 



NBpjuj^uj^B^m) [A(qu 2 )A(q'u 3 ) - A'(qu 2 )A'(q'u 3 )} 
(oj 2 —oJi — w 3 )(o; + u 2 - wi - tu 3 + irf) 



+^E/4(kqq') 



qq 




duidu 2 du>3 



N B (u 1 u 2 u 3 )B(k'tUi)B(qtu 2 )B(q'u 3 ) 
(oj 2 — ui — oj 3 )(uj + tu 2 - ui - u 3 + irf) 

+ A TE^(kq) JJ~ ' du,duj 2 

NpjuiUz) [A(k + q, u^Ajqu^ + A'(k + q, u^A'jqu^} 
(uo 1 - u 2 )(u + Ui - u 2 + irf) 



(8) 



where 



/l(kqq') 



</?i(q - k, k') - ^i(-q' - k, k') + -^ 2 (-q' - k, q) + ip 2 {-q' , q) 



5 



3 3 1 

-^(-q', q - k) + -¥> 2 (-q, q' + k) - -y 2 (k - q, q') 

-</? 2 (-q,q') , 

/ 2 (kq) = ^(_ q - k, k) - -^i(-q - k, 0) - -<^(q, 0)] , (9) 

1 11 
/ 3 (kqq') = <Pi(-q, q - q') ~ ~ k, k') - -y?i(q - k, k') + -<^ 2 (k - q, q') 

+ -<p 2 (-k - q', q) - -<£> 2 (-q', q - k) - -<£> 2 (-q, q' + k), 

/ 4 (kqq') = <f 3 (q, q' - q) - -<f 3 (-q, k + q') - ^(q', k - q) - ~</?s(q - k, k') 

+ 7^s(-k - q', q) + -p 3 (q - k, -q') - ~^ 3 (-k - q', k'), 

k' = k- q+ q', 7] = +0, 

^l(qq') = Yl I 6 ™ ~ Jj) Vp+q' = Mq+q' - 4t2 7q' ~ ^t'fjj, 

^2(qq) = \ 6 pi ~ Jj) f P J P+q' = ^q J q+q' - 4tJ7 q /, (10) 
p ' 

V?3(qq') = ^ Upq - J JpJp+cl' = <Vq+q' ~ 4 

N bf (ujiuj 2 uj z ) = [1 + n F (w 2 ) [1 - n F (u; 3 )] - n B (wi) [1 - " F (a; 2 )] n F (u 3 ), 

N F (u)iU) 2 ) = n F {ui) -n F (u) 2 ), (11) 
N b (uiu 2 ujz) = [l + n B (a;i)]n s (a; 2 ) [l + n B M] —n B (wi) [1 + n B (u 2 ))n B (uj 3 ). 

Here n B {oj) = {e w ^ T — l) 1 , n,ip(u;) = (e^/ 7 " + l) 1 with the temperature T, the Fourier 
transforms of the hopping and exchange constants are 

*k = 4t 7 k + 4i% J k = 4J 7k , (12) 

t' is the next nearest neighbor hopping constant and 7 k = cos{k x ) cos(fc y ), 

B(ku) = -7r" 1 ImD(kw), 

AQiu) = -7r- l ({a kl7 \al)), (13) 
A'(kw) = -n- 1 ((a^\a_ kl )) 

are the spectral functions of spin excitations and holes, ((ako-|o k(T )) and ((ctktl a -k4,)) are the 
retarded anticommutator normal and anomalous Green's functions. 

In the derivation of Eq. (jSJ) we have taken into account that due to the symmetry of 
Hamiltonian (pQ) 

((s^\s k 1 )) = 2D(ku J ). 
Notice that the polarization operator (JSJ) satisfies the following relations: 

Ren(kw) = -ReII(k, -u), Imn(kw) = Imn(k, -u) 
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and hence 

B(ku) = -B(k,-w). 

For to > and zero temperature the first term in the right-hand side of Eq. (jSJ) describes 
the contribution of the process, in which a spin excitation decays to yield a fermion pair 
and another spin excitation, the second term corresponds to the decay with the formation 
of three spin excitations, and the third term is connected with the transformation into a 
fermion pair directly and with assistance of a hole [see Eqs. ()A.14j) and flA. 17j) ] . 

In the below calculations, the following approximate expression for the spectral function 
of spin excitations is used in the integrands of Eq. (JSJ): 

BQxu) = 2J ^° 1 ^ 1 ~ 7k ) [5(u - u k ) - 6(u + u*)) . (14) 
w k 

This expression is derived from Eqs. (j3J) and §5§ if the damping of spin excitations is ne- 
glected to make tractable multiple integrations in Eq. (|H1). In Eq. ( fl4l) . the term a;Ren(ku;) 
in the denominator of Eq. (Tjf]) is also neglected, and the frequency of spin excitations is 
approximated by Eq. (E]), in which terms containing parameters F( are dropped for the 
considered small x. The presence of holes exerts main influence on the parameter 5.— In 
the case of the short-range antiferromagnetic order, which is settled under the influence of 
thermal fluctuations 43 and/or a finite hole concentration, the spin-excitation dispersion is 
gapped at Q (see Fig. [5] below). The parameter 5 defines the value of this gap, which is 
directly connected with the spin correlation length. Values of the parameters in Eq. (JTJ 



were taken from the self-consistent calculations of Ref. |44. 

In the superconducting state, the hole spectral functions for wave vectors near the Fermi 
surface are approximated by 



ah ^ + e k ) Z(£ k — e k ) . Z 

A{ku) = — 5{ U -h)+ — Siu + ^ + ^-^Siu) 

ZA. 

A'(ku) = —± [6{u - &) - S(lj + &)] . 



(15) 



where 



1 if Eq - 2t < u < -A, 
S(u)={ or A < u < e Q + 6t, (16) 
in other cases, 



e k and ^k = \Ak (^ k ) 2 are ^ e ^°^ e dispersions in the normal and superconducting states, 
respectively, with the ci-wave superconducting gap function A k = A[cos(k x ) — cos(fc 3/ )]/2. 

The function S(u) in A(ku>), Eqs. (|T5|) and ( TT6|) . describes the strong HSE continuum, 
which is typical for the hole spectral function in the t-J model. The continuum spans the 
frequency range approximately equal to the width of the initial band, in the present case 
81 For moderate doping the lower edge of the continuum lies approximately 2t below the 
bottom of the hole band e^,— which is located at the momentum Q in the used hole picture. 
In accord with the definition of the spectral function the spectral weights of its coherent and 
incoherent parts are connected by the relation 

z + z'= l -±^. 
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Below we set J = 0.2t. This ratio of the parameters belongs to the range derived for hole- 
doped cuprates.— For these parameters and hole concentrations 0.07 < x < 0.15 the ratio 
between the spectral weights is Z/Z' ~ 0.1.— 

We have approximated S(u) by constant values within certain frequency ranges to sim- 
plify integrations in Eq. (jSJ). In Eq. (TT5|) . the intensity of the HSE continuum is set equal to 
zero in the interval —A < u < A. For the considered c?-wave order parameter this intensity 
may be nonzero in the ranges —A < u> < — |Ak| and |Ak| < u> < —A. However, since the 
hole damping decreases as the Fermi level is approached, it can be supposed that the contin- 
uum intensity in these ranges is much smaller than far apart from u — 0. Notice, however, 
that its perceptible intensity in these intervals may be of critical importance for the charac- 
ter of the magnetic response below the resonance frequency. In the calculations we increased 
the continuum intensity in these regions up to the value it has far from the Fermi level. The 
response became commensurate due to the growth of |ImZ)| at the antiferromagnetic wave 
vector. 

For the considered hole dispersions corresponding to moderate x and for frequencies 
u < 0.3t = 150 meV the main contribution to the magnetic susceptibility from coherent 
peaks of the spectral functions (fl5|) is made by the last term in the right-hand side of Eq. ([8]) 
(here and below in going to energy units we proceed from the estimate^ 1 ^ t = 500 meV). 
As mentioned above, this term is connected with the process of the spin-excitation decay 
into fermion pairs. The substitution of Eq. (fT5l) into ([8]) gives an addend with the multiplier 
5(u> — £k+q — £q) i n the imaginary part of this term. The approximation used for the 
representation of this 5-function has an important impact on the character of the magnetic 
response below the resonance frequency. In fact, this approximation is connected with the 
shape of the coherent peaks in the hole spectral function. If this 5-function is approximated 
by a Lorentzian with a constant width T, the low- frequency incommensurate response arises 
only for a Fermi surface with nesting, as it was noted for the analogous term in the itinerant 
electron theory— However, not all dispersions derived from the photoemission of p-type 
cuprates possess nesting. Nevertheless, the incommensurate magnetic response is observed 
also in crystals with such dispersions. As indicated above, the problem is connected with 
far tails of the Lorentzian. They strengthen Imx at the antiferromagnetic wave vector. This 
hides incommensurate peaks, which are less intensive in the absence of nesting. However, 
the decay of tails of the spectral peaks has to proceed faster than in Lorentzian, since the 
damping of the respective states grows with distance from the Fermi level. Besides, the 
tails have to be cut off in the superconducting gap. To take these facts into account, we 
approximate the function 5(u — £k+ q — £q) by the n-shaped step function 



With this approximation, the low-frequency incommensurate response was obtained for sev- 
eral hole dispersions used earlier for p-type cuprates,— including those derived from 
photoemission data. 

For cu > and T = the result of the substitution of spectral functions and (|T5|) 
into Eq. ([H]) is given in Eq. flA. 17[) . which was used in the present calculations. 




(17) 
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FIG. 1: The Fermi surface of the dispersion fj 18|) . 



III. RESULTS AND DISCUSSION 



In this section, figures were obtained with the following dispersions for the coherent part 
of the hole spectrum: 



£k 
t 



cos(k x ) + cos(ky) — 0.3 x cos(k x ) cos(k y ) — 0.2, 



(18) 



-0.1305 + 0.5951 x 



cos(k x ) + cos(ky) — 0.1636 x cos(k x ) cos(k y 



+0.0519 x - 



+ 0.1117 x - 



cos(2A; a ;) + cos(2k y 
+ cos(k x ) cos(2k y ) — 0.0510 x cos(2/c ;c ) cos(2A;j / ) 



cos(2/c x ) cos(ky) 



(19) 



Dispersions similar to that given by Eq. ffTHl) a re frequently used for the description of hole 
bands in p-type cuprates (see, e.g., Refs. 20| 48 50). Its Fermi surface is shown in Fig. ED The 



second dispersion (|T9|) was borrowed from Ref. [21 
data in Bi 2 Sr 2 CaCu208. Notice that if band (EE 



It was derived from the photoemission 



were 
The 



and similar bands in Refs. |20|J48|J50 
used for underdoped cuprates, x < 0.12, the dispersion (|T9|) corresponds to x ~ 0.17.— 
coefficients in Eq. ( IT9|) are in electronvolts. In this equation, signs are opposite to those in 
Ref. |2l|, since we use the hole picture. The general shape of the Fermi surface for dispersion 
(1T9"|) is similar to that shown in Fig. [T] 

The other parameters used in the below calculations are the following: J/t = 0.2, Z/Z' = 
0.1, T = 0.002t, 5 = 0.25, A = 0.05t, C x = -0.2, and a = 1.7. The last two parameters 
are taken from the self-consistent calculations for low hole concentrations^ and are close to 
the results for the undoped case.— The parameter 5 is somewhat increased in comparison 
with the self-consistent result of Ref. |44j to compensate the correction introduced by the term 
a;ReII(ku;) in the present calculations. We chose for 5 the value, which provides the spin gap 

Wq = [oJq + WQReII(Qu;Q)] ^ 2 ps 40 meV with dispersion ( TT8|) . Gaps of such a magnitude 
were obtained in the self-consistent calculations for moderate hole concentrations— This 
value is also close to the frequency of the resonance peak observed in some cuprates.- 1 ^ 
Results of calculations with Eqs. (jl]) and flA.17j) are given in Fig. [2] for dispersion f|T8|) and 
in Fig. [3] for dispersion (EL9l . 
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FIG. 2: The imaginary part of the spin Green's function for the hole dispersion (|18|) : uj = 0.04i = 
20 meV (a), 0.077i = 38.5 meV (b), and 0.12t = 60 meV (c). Other parameters are given in the 
text. 



As seen in the figures, for both dispersions the magnetic response is incommensurate for 
small frequencies (Figs. [2(a) and [3(a)), and the maxima are located on axes of the Brillouin 
zone. With increasing frequency the response becomes commensurate (Figs. E(b) and [3(b)) 
and reaches the maximal intensity. For even larger frequencies, the response is again incom- 
mensurate (Figs. [2(c) and [3(c)). However, this time it looks more isotropic than for small 
frequencies. Similar behavior of the magnetic susceptibility are observed experimentally in 
p-type cuprates (see, e.g., Ref. |l2|). For frequencies u > 0.26t = 130 meV for dispersion f [T8|) 
and u > 0.2t = 100 meV for dispersion f[T9|) with the used set of parameters the momentum 
dependence of the susceptibility transforms to a broad commensurate maximum. For both 
dispersions for frequencies u < 0.02t = 10 meV the susceptibility is peaked at the diagonals 
of the Brillouin zone. 

Notice that in Ref. |2l| no incommensurability was found for small frequencies with disper- 
sion (fT9l) . in contrast to the result shown in Fig. [3(a). As mentioned above, this discrepancy 
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FIG. 3: The contour plots of — t!mD(ku}) for the hole dispersion f)19[) : oj = 0.06i = 30 meV (a), 
0.0747i = 37.4 meV (b), and O.li = 50 meV (c). Other parameters are given in the text. 



is connected with different representations of the ^-functions ensuring the energy conserva- 
tion in the convolutions of coherent parts of the hole spectral functions. With representation 
( I17p results similar to those shown in Figs. [2] and [3] were also obtained for some other bands 
used earlier for the description of p-type cuprates. Susceptibilities obtained with different 
dispersions vary somewhat in the location and intensity of maxima and in some other de- 
tails. Such discrepancies can be expected, since the dispersions were obtained for crystals of 
different compositions and with dissimilar hole concentrations. For the bonding band from 
Ref. |49| we found low-frequency maxima on diagonals of the Brillouin zone, while for other 
considered bands they are on the axes, at least for not too small frequencies. 

For band (Il8p the dispersion of susceptibility maxima, which are located on the X axis 
of the Brillouin zone, is shown in Fig. HJ It is similar to the hourglass dispersion observed 
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FIG. 4: The dispersion of the susceptibility maxima along the axis of the Brillouin zone for the 
parameters of Fig. [2j The momentum is in the reciprocal lattice units. 

experimentally.— Notice that in Fig. H] the wave vector is in the reciprocal lattice units, 
2ii/a = 1 r.l.u. As seen in the figure, for dispersion ([IB]) the maxima are located at k x = 
at low frequencies, which is typical for the hole concentration x = -.— 

The formation mechanisms of the incommensurability below and above the resonance 
frequency are different. Green's function (j4|) has the resonance denominator, which is re- 
sponsible for the incommensurate response above u r . If we drop temporarily the real part 
of the polarization operator in the denominator of Eq. (j4j) the resonance condition reads 

u = u k . (20) 

As follows from Eq. fl6]), near Q the spin-excitation dispersion is given by the expression 



« JwJ, + c 2 (Q - k) 2 , c=y/8a\C 1 \J. (21) 

This dispersion is shown in Fig. |5] As mentioned above, the dispersion has a gap at Q 
and the magnitude of the gap uq is connected with the correlation length of the short- 
range antiferromagnetic order.— « For u > uq the graphic solution of Eq. ( 120]) is shown 
in Fig. |5j Thus, if the imaginary part of the polarization operator is not too large, the 
susceptibility is peaked at momenta indicated by arrows in the figure, i.e. on the circle 
surrounding Q, 37,38 as seen in Figs. Hfc) and[3](c). This is qualitatively the mechanism of 
the formation of the incommensurate response above the hourglass neck. The inclusion of 
Reil(ka;) leads to a renormalization of ojq and dispersion ( 1211) immediately above the spin 
gap. For larger frequencies the real part of the denominator in Eq. becomes nonzero. 
Besides, the damping ImIT(ka;) grows rapidly with u. As a consequence the high-frequency 
incommensurate maxima are less pronounced or a broad commensurate peak is formed. 
Such shapes of the susceptibility were also observed experimentally.— 

Below the resonance frequency the momentum dependence of the imaginary part of the 
susceptibility is determined by Imll(kx>) in the numerator of the expression 

[oo 2 - wRen(kw) - ulf + [wimnxMf 



ImD(kc) = ' tk " lllil V ; 2- (22) 
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k=Q 

FIG. 5: The dispersion of spin excitations © near the antiferromagnetic wave vector Q. Wave 
vectors, for which the resonance condition (|20p is satisfied, are shown by arrows. 

\ 

/ 

\ 

y 

-6 -4 -2 2 4 6 

K 

FIG. 6: The contour plot of the dispersion £jj = + A^, where is given by Eq. (fT8|) and 
A = 0.05i. Contours correspond to the energy 0.03i. The squares shown by thin solid and dotted 
lines are the first and magnetic Brillouin zones. The blue dashed lines and the red dotted line 
connect fermion pairs contributing to the susceptibility on the axis and diagonal of the momentum 
space, respectively. 




In this frequency range, the last sum in Eq. (J8]) [which corresponds to the last sum in 
Eq. (1A.17P ] makes the main contribution to the polarization operator. This sum is connected 
with the decay of a magnetic excitation into two fermions. The incommensurate peaks in the 
susceptibility are related to the contribution of coherent parts of the hole spectral functions 
into this sum [the term with the multiplier P(u — £k+q — £q) m the last sum of Eq. (IA.17j) ]. 
The fermions with energies £k ~ w /2 play a key part here. The location of these states in 
the momentum space is shown in Fig. [61 In this figure, crescent pockets are sections of the 
fermion dispersion £k = \A-k + a ^ the energy 0.03t. For small energies this dispersion 
looks like flattened cones originating at the node points. The blue dashed and red dotted 
lines connect pairs of fermion states, which give the main contribution to the magnetic 
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susceptibility on an axis and diagonal of the momentum space, respectively. The location of 
these states is restricted to the crescent pockets, and the difference in their wave vectors is 
equal to the spin excitation momentum. As a consequence the low-frequency susceptibility 
is sharply peaked at certain regions of the Brillouin zone, as seen in Figs. |2fa) and[3^a). The 
line segments shown by dashed and dotted lines in Fig. E] differ from the antiferromagnetic 
wave vector and, therefore, the response is incommensurate (in this figure, opposite sides 
of the magnetic Brillouin zone are displaced by the vector Q from each other). This is the 
mechanism of the formation of the incommensurate response below the resonance frequency. 
Notice that it implies the fulfillment of the energy conservation: oo = £ q + £k+q [see the 
mentioned term in Eq. flA.17j) ]. In the processes, in which distant tails of the coherent 
peaks contribute to the polarization operator, this energy conservation is violated. If this 
contribution dominates the low-frequency response will be commensurate. 

For frequencies immediately below oo T a number of fermion pairs, which contribute to the 
susceptibility peak on an axis of the Brillouin zone, is larger than for the maximum on a 
diagonal. This conclusion can be drawn from Fig. [5]- there are two groups of fermion pairs, 
which make the main contribution in the former case, and only one group in the latter. 
This explains the maxima of — lmD(ku>) on the axes in this frequencies range. The same 
location of susceptibility maxima is observed in experiment.-^ - - However, with decreasing 
co the size of the crescent pockets is reduced and the intensity on axes goes down faster 
than on diagonals in our calculations. Finally, at very small frequencies the susceptibility 
becomes peaked at the diagonals of the Brillouin zone (for parameters of Figs. [2] and [3] it 
happens at uo = 0.02t). For even smaller frequencies the intensity on the axes disappears at 
all. The analogous behavior was observed in the itinerant-electron approach.— This result 
seems contrary to experimental data, where the susceptibility is peaked on the axes up to 
to = 0.— ~— A possible way to resolve this contradiction is discussed in the next section. 

With increasing the hole concentration regions of the Brillouin zone, which are filled by 
holes, grow.— These regions are revealed by the crescent pockets around the zone corners 
in Fig. [6j As a consequence dashed and dotted line segments in the figure move more and 
more away from Q. Thus, the incommensurability parameter k grows with x. This growth 
is also promoted by the concentration dependence of the spin-gap frequency cjq. As follows 
from Eq. f[2"2"j) . for small frequencies 

ImD(ku;) « w/i k o;^ 4 Imn(ku;). 

In this equation, the multiplier ui^ 4 is peaked at Q [see Fig. [5] and Eq. (1211) ] and therefore 
it impedes the appearance of the incommensurate peaks. However, the gap cuq grows with 
x,— effectively slackening the momentum dependence of cu^ 4 and promoting the separation 
of the susceptibility maxima from the antiferromagnetic wave vector. The indicated increase 
of the incommensurability parameter with x is observed experimentally. 3 

As seen in Figs. [2] and [3l for the used parameters in the region of the hourglass neck 
the susceptibility peaks are much more intensive and narrower than the maxima for smaller 
and larger frequencies. This peculiarity of the spin-excitation spectrum is called the reso- 
nance peak.— The formation mechanism of the peak is explained in Fig. [3(b), in which the 
real and imaginary parts of the denominator of the spin Green's function (jlj) are shown. 
The considered crystals have two energy parameters, which may have comparable values. 
The first of these parameters is the spin-gap frequency uj'q, for which the real part of the 
mentioned denominator vanishes at k = Q. The second parameter is the lower edge of the 
two-fermion continuum, cu e (Q) = min(£ q + £q+q) at the antiferromagnetic momentum. Due 
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FIG. 7: (a) The resonance peak in the spin susceptibility, (b) The real (the blue solid line) and 
imaginary (the red dashed line) parts of the denominator of Eq. Q . Parameters are the same as 
in Fig. [2j except that 6 = 0.16. 



to a small but finite damping of fermion states this edge appears as a fuzzy step near 0.08t 
in the imaginary part of the denominator in Fig. [7(b). This step has the same origin as in 
the theory of itinerant carriers,— i2£i2I where for a vanishingly small damping it transforms 
to a discontinuous jump leading to a logarithmic singularity in the real part of the nonin- 
teracting susceptibility. In this theory, the resonance peak is related to this singularity. As 
seen in Fig. UJb), the small damping T = 0.002t = 1 meV transforms the singularity into 
the weak maximum. For the parameters of Fig. [7J as well as for Figs. [2] and El < ^e(Q), 
the real part of the denominator of Eq. (J4j) vanishes in the region of small damping below 
w e (Q), which produces the narrow and intensive resonance peak.— Apparently in this case 
bj T coincides with the spin-gap frequency u'q. It grows with x, as observed experimentally 
in YBa2Cu307_^, 9 since with increasing the hole concentration the spin correlation length 
decreases. Notice that the weak damping below w e (Q) is mainly connected with the contri- 
bution of the HSE continuum [the second and third terms in the last sum of Eq. (IA. 1T[) . 
As the wave vector recedes from Q the spin-excitation frequency crosses the edge of the 
electron- hole continuum w e (k), which separates regions of weak and strong damping. The 
incommensurate peaks become broader and less intensive and finally the magnetic response 
becomes commensurate, as mentioned above. The considered case resembles the situation 
in YBa 2 Cu307_.y and Bi 2 Sr 2 CaCu208+<5.-~- 

Another situation occurs if u'q > co e (Q). Due to a larger damping the susceptibility 
maximum at Q is broader and less intensive than in the previous case. In particular, this 
situation takes place for A = 0.03t if the other parameters are the same as for Fig. [2J Due 
to broader maxima the hourglass dispersion looks somewhat different than that shown in 
Fig. H] - the neck with commensurate susceptibility maxima occupies a substantial range 
of frequencies (for the mentioned parameters from 0.05t to 0.12t). The broad maximum 
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with a moderate peak intensity at the antiferromagnetic wave vector resembles the picture 
observed in lanthanum cuprates.— 1 ^ 

IV. CONCLUSION 

The above discussion was referred to the superconducting state. To consider the pseudo- 
gap phase, we have to drop the contribution of the anomalous hole Green's functions and 
to introduce a finite damping in the pseudogap region \cj\ < A. As mentioned above, the 
damping blurs the incommensurate peaks or even leads to a commensurate response. Such 
changes in the susceptibility are observed experimentally above T c .— Additionally it should 
be taken into account that in the pseudogap phase, the depression of the hole spectral in- 
tensity near the Fermi level is not described by the (i-wave function. In this state, near the 
diagonals of the Brillouin zone gapless regions of the Fermi surface - arcs - are located.— 1 ^ 
As mentioned above, in the superconducting phase with decreasing frequency the intensity 
of the susceptibility maxima on the axes goes down more rapidly than on the diagonals. As a 
result for small frequencies the intensity is peaked at diagonals. This fact is connected with 
shrinking crescent pockets of fermion states, which contribute to the susceptibility maxima, 
when oj —> 0. Obviously, with the appearance of arcs in the pseudogap phase the situation 
is changed: with decreasing frequency the pockets do not shrink to points on the diagonals 
but rather acquire the shape of the arcs. Consequently, it can be expected that the sus- 
ceptibility maxima will be located on the axes up to zero frequency. The incommensurate 
elastic response with maxima peaked on the axes is observed in lanthanum and yttrium 
cuprates.— ~— Based on the sharp difference in the spin correlation lengths below and above 
the energy E g ~ 4 meV it was supposed^ that below T c the crystal Lai.875Sro.i2sCu04 is 
phase-separated into superconducting and non-superconducting regions. From the above dis- 
cussion one can expect that the incommensurate elastic response is due to the latter regions, 
which are in the pseudogap state, and both regions contribute to the inelastic susceptibility 
with the hourglass dispersion for u > E g . This picture allows one to suggest a mechanism 
for the strengthening of the incommensurate elastic response by magnetic fields H < H c2 , 
observecr^ 1 ^ 1 ^ in lanthanum and yttrium cuprates. The growing field increases the part of 
the crystal occupied by vertex cores, which are in the pseudogap state and contribute to the 
response. 

In summary, we calculated the magnetic susceptibility of the two-dimensional t-J model. 
To consider the case of strong electron correlations the Mori projection operator technique 
was used. All processes contributing to the spin polarization operator - the decay into a 
fermion pair, the decay into the pair assisted by a hole or with a creation of a spin excitation, 
and the decay into three spin excitations - were taken into account. In these calculations, 
both the coherent and the strong incoherent parts of the hole spectral function were allowed 
for. For the coherent part several dispersions found in the literature were analyzed, including 
those designed to fit photoemission data. For the superconducting state with the c?-wave 
order parameter we found that at certain conditions all considered dispersions produce the 
incommensurate magnetic response below the resonance frequency. For this incommensura- 
bility the nesting of the Fermi surface is not required. This result proposes the explanation 
for the fact that the low-frequency incommensurability is observed in a wide range of hole 
concentrations and in different families of p-type cuprates. The mentioned conditions re- 
duce to the fast decay of the tails of the hole coherent peaks and to a small intensity of 
the hole-spin-excitation continuum near the Fermi level. Both conditions are related to the 
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damping of hole states, which decreases on approaching the Fermi level. The used equations 
allowed us to obtain the hourglass dispersion of the susceptibility maxima. Their locations 
in the momentum space are close to those observed experimentally. The upper part of 
the hourglass dispersion is related to the dispersion of spin excitations. The part below 
the resonance frequency is connected with the incommensurate maxima of spin-excitation 
damping, which is mainly caused by the decay into a fermion pair in this frequency range. 
Immediately below the resonance frequency the susceptibility is sharply peaked at the axes 
of the Brillouin zone due to a larger number of final states in comparison with the diagonals. 
The intensity distribution in the upper part of the hourglass dispersion is more isotropic. 
In this picture, the growth of the low-frequency incommensurability parameter with doping 
is related to the increasing area occupied by holes in the Brillouin zone and to the rising 
gap in the spin-excitation dispersion at the antiferromagnetic momentum. The resonance 
peak is narrow and intensive if the gap frequency falls on the region of small damping below 
the edge of the two-fermion continuum at the antiferromagnetic momentum. In this case, 
the resonance frequency is close to the gap frequency. If the gap frequency exceeds the 
edge of the two-fermion continuum, a broader and less intensive commensurate maximum 
is produced. 



Appendix 

In this Appendix equations of Sec. II are derived. In the Mori projection operator tech- 
nique D(ku) is calculated from the Kubo relaxation function 

/*oo /»oo 

((si\s\))= dte** ^([4(*Vi k ]>. (A.l) 
Jo Jt 

The two functions are connected by the relation 

u{(sl\st k )) = (s* k ,s\)+D(ku), (A.2) 

where 

/■oo 

(st,s\)=i dt([s z k (t),s\]). (A.3) 
Jo 

In this approach, Green's and relaxation functions are represented by continued fractions. 
The elements of these fractions can be calculated using the recurrence relation, 5 ''' which 

resembles that of the Lanczos algorithm.— As applied to the relaxation function (lA.ip . the 
continued fraction and the recurrence relation read 

((4K k )) = , (A.4) 



u-Ei- V^ku) 
[A n , H] = K-iA>-i + E n A n + A n+1 , E n = ([A n , H],A n ) (A n , A^y 1 , 



V n = (A.+1, K, A{) 1 , V-! = 0, A = S 



(A.5) 



k , 71 = 0,1,2... 



where the inner product of operators (A, B*) is determined by Eq. f]A.3|) . The operators A n 
obtained from the recurrence relation (1A.5[) are orthogonal, 

A n ,Al) oc<W- (A.6) 
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The residual term R(ku) in Eq. (1A.4I) is given by the relations 

-l 



R(kt) 



A 2t ,4) [A 2} A 



it 

2) ^2 



i?(kw) 



dte wt R(M), 



where the time dependence of the operator A 2t is determined by the equation 

d 



1 - P )(l - ^2,t=0 = v4 2 , 



(A.7) 



(A.8) 



with the projection operators P n defined by the relation P n Q = (Q, Ajj) (A n , A^) 1 A n . 

From Eq. (|A.5[) and the definition of the inner product flA.3[) it can be seen that in the 
continued fraction (IA.4I) 

E = Ei = 0. (A.9) 
Using these results in Eqs. flA.2j) and (1A.4I) we obtain Eq. (T4J). In these equations, 



& k = (a, 4) = <K,-^-J) 



with zsi k = [si k , i?] and 

n(kw) = Ki2(kw), ^ = 1/ - (A.10) 

To calculate Green's function (j3J) we have to estimate quantities (lA.lOj) . In accordance 

with Eqs. (I A. 51) and (|A.9j) . the operator A2, which enters into the polarization operator 
II(ku;), is equal to 

A 2 = i 2 sl-ulsl (A.ll) 
The quantity o; k is given by the relation 



s k5 s -k) 



-1 



(A.12) 



The second term in the right-hand side of Eq. (1A.11I) ensures the orthogonality ( 1A.6I) of the 
operators A = s k and A 2 . Notice that the orthogonality of the operator 



(A.13) 



M - ist - TTTff y2 y2 ( tk ' ~~ ^k'-k) a^ak'-k^ - (A> - Jv-u) s k' ls k'-k 

2 V N V 1 rr J 



to the operators A and A 2 is provided by conditions ( 1A.9j) . Thus, the calculation of A 2 is 
reduced to the separation of the quantity i 2 s^ into a part proportional to s k and a part or- 
thogonal to this latter operator, in accord with the Mori orthogonalization procedure. 36 This 
separation cannot be performed exactly, since the quantity (s k , sl k ) in Eq. f ]A.12j) cannot be 
directly calculated. To carry out the separation approximately we use the decoupling, which 
is similar to that applied for the Heisenberg model.— For the t-J model the analogous 



approach with some additional approximations was used in Refs. |39|J44|J60 
The operator £ 2 s k reads 



Z 2 5 k 



^ A( k( l<l> s -ak' a l-aa q v + -7= h( kc fi aa l+<i,<T 

qq'cr V qcr 



"qcr 



+ AT /3(kqq')4' a qa a qv + U (kqq) 4' S+ 1 s q , x 

qq'cr qq' 
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+ T7 (*q+q' — ^q-k+q') S k-q'^q'' 



N 

qq' 



+^ E J o ( " J q-k+q') 4- q '*?, ( A - 14 ) 



N 

qq 



where the functions /<, t q , J q are given by Eqs. © and QEJ, X°° = A^" 1 ^ e"^ 1 ^ 00 , 
x tt = TV" 1 ^e-^ and = 

The above equation for i 2 s^ needs in some explanation. In this equation, terms containing 
hole operators can be expressed in two equivalent forms. In the site representation these 
forms read 

(X?° + X?°) a\ la a v , a (A. 15) 

and 

Q + l X i° + <™i) 4.°!"- (A.16) 
Their equivalence follows from the completeness condition of the site states, 

However, after the decoupling used for calculating II(ka;) this equivalence is lost. If 



Eq. (IA.15I) is used, as it was done in Ref. |39|, the decoupling of the operators X[ a leads 



to a small contribution of the decay into a fermion pair, which is assisted by charge fluctu- 
ations. In this case, the contribution of the direct decay of a spin excitation into a fermion 
pair is gone. As follows from Eq. (1A.16I) . this process exists. Among terms in the po- 
larization operator this is the simplest process, which makes the main contribution in the 
most interesting range of frequencies u < 0.3t = 150 meV. As mentioned above, the in- 
commensurability below the resonance frequency stems from this term. Besides, neglecting 
the decay into a fermion pair is in contradiction with the spin-wave and itinerant-electron 
approximations for the susceptibility, in which only this process is usually considered. Thus, 
representation (1A.15I) is not applicable to the considered case. Equation (1A.14|) is derived 



using (1A.16j) . This expression allows one to take proper account of the decay into a fermion 
pair, gives simpler decoupling and does not introduce additional Green's functions. 

In all terms of Eq. flA.14j) operators belong to different lattice sites. The part of the 
operator 2 2 sj<, which is proportional to s^_, is contained in the last four terms in the right- 
hand side of this equation. This part is obtained by substituting operator multipliers of s£ 
in these terms by their mean values. From these mean values the quantity u^, Eq. ([H]), is 
formed. From Eqs. §5§, (JSJ) and (1A.12j) the quantity fs k , sf_ k ) can be derived, 

( s k> s -k) 

In the absence of holes this equation, (jBJ) and (JZj) reduce to the formulae obtained earlier for 
the Heisenberg antiferro magnet.— & 

The operator is obtained from i 2 s^ by subtracting terms entering into u^s^.. The 
contribution of the two last terms in Eq. (IA.14j) into can be neglected, since charge 



J|Ci 


-tF x 


AJ(Ja 


Ci\ 


-tF 1 )(5 + l + -y k ) 



fluctuations are small in the t-J model for small hole concentrations. Thus, A 2 contains two 
first terms of Eq. (IA.14I) together with the third and fourth terms of this equation, in which 
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a qo- a q'o- an d s q ls q' are replaced bya* a q v — (a qcr a q / CT ) and Sq 1 ^/ 1 — (Sq 1 ^ 1 ), respectively. 
Substituting A 2 into Eqs. (1A.7[) . (lA.lOp . approximating the operator A 2 t by A 2 (t) and using 
the decoupling, we obtain Eq. flS]). 

Inserting Eqs. ffT4"|) and f|T5|) in (jSJ), for u; > and T = we get the equation used in the 
above calculations, 



1 f4vr|Ci|J 



ImIT(M = -^{^^E[ 2 A 2 ( k ^') + /3(kqq' 

qq' 



ai 1 - 7k' 



4 4q4q' 



4£ q (4t - A) 



■6(uj - £ q / - uv - A) + 



8t- 2A 



4^(4* -A) 

32tt(J|C 1 |) 3 ^ (1 - 7k ,)(l - 7q )(l - TqQ 

^ 4 o; k /o; q a; q / 



U(oj — uj^i 



N 2 



qq 



X?(Wk' + W q / + W q - W) 

Z 2 [(e k+ q - £k+q)(£q + £q) " A k+q A q ] 



+f £/ 2 2 (kq) 

q 

xP(w-£ k+q -£ q ) + 



4£k+q£q 



4£ q (4t - A) 
where k' = k — q + q' and 



4£ k+q (4t - A) 



Z' 



8t-2A 



U(u) 



(A.17) 



I7(w) = 

in the considered range of frequencies 



if a < 2A, 
u - 2A if 2A < u < 2t - e Q + A 
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